function [ psiK ] = createPsiK( kx0,dkx0,x0,L,nm )
    eps = 0.000;
    n2 = nm / 2;
    
    A = [-n2:n2-1];
    A = A*(2*pi/L);
    
    psiK = zeros(nm,1);
    for k =1:nm
        ck=  getCk( A(k),kx0,x0,dkx0);
        if(abs(ck)> eps)
            psiK(k) = ck;
        end
    end
       
end

function [ ck ] = getCk( kx, kx0 , x, dkx)
    dkx = dkx * 2 * pi;
    coeff =power( prod(dkx),1/4 );
    real = ((kx - kx0).^2)./(dkx * 4);
    ck =  (exp(-sum(real))/ coeff) * exp( -1i *kx*x);
end
    